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We calculate the equation of state of pure neutron matter, comparing the G-matrix calculation 
with the in-medium T-matrix result. At low densities, we obtain similar energies per nucleon, 
however some differences appear at higher densities. We use the self-consistent spectral functions 
from the T-matrix approach to calculate the ^5*0 superfluid gap including self-energy effects. We 
find a reduction of the superfluid gap by 30%. 
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In the interior of neutron stars, neutron rich nuclear 
matter is present in a range of densities up to several 
times the normal nuclear density. Two properties of the 
neutron matter are important for the modeling of neu- 
tron stars. The equation of state, which serves as an 
input for the calculation of the mass of the star, and 
the superfluid parameters of the dense matter, relevant 
for the cooling and the rotation of stars. Dense neutron 
matter cannot be studied directly in laboratory; there- 
fore, it is important to have a reliable description of 
this many-body state from theoretical calculations. Very 
advanced methods based on the Brueckner-Hartree-Fock 
(BHF) approximation were applied for the calculation of 
neutron matter properties . Independently, approaches 
using the variational chain summation to tackle the nu- 
clear many-body problem |§, || have been developed to 
a high level of accuracy. For symmetric nuclear matter, 
these calculations have shown the need to introduce a 
three-body force to reproduce the phenomenological sat- 
uration point. Relativistic corrections have also been es- 
timated 1^ . Variational and BHF methods allow to study 
the equation of state of neutron matter up to densities 
several times the normal nuclear matter density. Micro- 
scopic calculations of different pairing gaps in neutron 
matter are also available ^, ^ . 

An alternative many-body scheme for the calculation 
of the properties of nuclear matter is the self-consistent 
in-medium T matrix |[| 0, |ll) 



T^V + VGGT 



(1) 



where the Green's functions G^^ — uj~p'^/2m — I] in 
the ladder propagator are dressed by the self-consistent 
self-energy 



iE = Tr[TAG] 



(2) 



The T-matrix approach is a conserving approximation 
1^ and gives consistent single-particle properties. The 
details of the equations and of the solution in the real 
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time formalism can be found in ^, Since the self- 
energy has a non-zero imaginary part, we obtain in the T- 
matrix approximation a nontrivial single-particle spectral 
function. 

At low densities, the T-matrix ladder diagrams are 
dominated by the particle-particle propagators and 
therefore should reduce to the G-matrix result (if the 
influence of the nucleon dressing is small at low densi- 
ties). However, at higher densities the ^-derivable T- 
matrix approximation takes into account some diagrams 
of higher order in the hole-line expansion. It is then 
instructive to compare the results of the T-matrix cal- 
culation with the corresponding BHF calculation for a 
realistic interaction. 

We perform self-consistent iterations of the T-matrix 
equations in pure neutron matter, using the separable 
Paris interaction fl^, |l5| and the same methods as in 
Ref. |l^. In Fig! ^ the resulting energy per nucleon 
is plotted as a function of the Fermi momentum. It is 
compared to the BHF calculation with the continuous 
choice of the single particle energy, using the same inter- 
action and numerical grids. An excellent consistency of 
the two methods is found at low densities, as expected. 
This provides a powerful check of the methods used for 
the solution of integral equations, which for the T-matrix 
case involve refined numerical algorithms. At higher den- 
sities the neutron matter in the T-matrix calculation is 
less bound, leading to a considerably harder equation of 
state. It is known that the BHF results with the contin- 
uous choice are not modified considerably by 3 hole-lines 
corrections The T-matrix approximation is a 

different summation of some of the 3 and more hole-lines 
diagrams. We obtain the first quantitative estimate of 
the difference between the self-consistent T-matrix ap- 
proximation and the BHF approach. Since the latter one 
is believed to be close to the converged result, we find 
an estimate of the accuracy of the T-matrix calculation. 
Below kp — 1.3fm^^ the difference in the energy per par- 
ticle from the two calculations is less than 2MeV. The 
discrepancy at higher densities could be attributed to 
the ring diagrams contribution. Also at higher densities 
meaningful results can be obtained only after inclusion 
of higher partial waves and three-body forces effects. In 
Fig. results of modern variational calculations |^ and 
of a BHF calculation including contributions from 3 hole- 
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FIG. 1: The energy per nucleon in pure neutron matter as 
a function of the Fermi momentum for the T-matrix (dotted 
line) and the BHF (sohd line) calculations using the Paris po- 
tential. The dotted line represents the results of a variational 
chain summation and the dashed-dotted line the results 
of a BHF calculation with 3 hole-lines contribution jl^ , both 
for a two-body Argonne V18 interaction. 



lines diagrams |Tj are also shown. Both methods are us- 
ing a two-body Argonne V18 potential and give similar 
results. As shown in Ref. ||T^ the continuous BHF results 
are almost indistinguishable from the full BHF calcula- 
tion with 3 hole-lines contributions included. Therefore, 
the difference between our BHF result (solid line in Fig. 
|l|) and the variational or next order BHF results must be 
attributed to a different interaction used, and to the lim- 
ited number of partial waves taken. The T-matrix energy 
per particle is larger than obtained in other approaches. 
Most importantly it is significantly larger than the con- 
tinuous BHF calculation with the same interaction, ft 
indicates that the T-matrix calculation does not give the 
correct energy per particle at higher densities. 

For the analysis of superfluidity, we use the spectral 
function obtained in the T-matrix calculation. As noted 
above, the T-matrix approach gives reliable results for 
the single-particle properties in the medium. The effec- 
tive mass in our calculation decreases from 1.07 to 1.02 
nucleon masses with increasing density, which is only 
about 0.05 more than the corresponding BHF effective 
mass. 

The superfluid gap equation with full spectral func- 



i2ny 



V{p, fc)A(fc) 



(3) 



where Ag{p,Lu) denotes the spectral function of the nu- 
cleon, including the diagonal self-energy T,{p,uj) (ob- 
tained in the T-matrix approximation) and the off- 
diagonal self-energy A{p) (obtained from Eq. (||) itself); 
A{p, oj) is the spectral function of the nucleon dressed 
with the diagonal self-energy only. Self-energy effects in 
the gap equation were studied in a number of papers 
[|l8|, |0[ |2|]. It was found that generally the super- 
fluid energy gap is reduced by the in-medium dressing of 
nucleons. The mechanism of this reduction can be under- 
stood easily from the effective quasiparticle gap equation 

MM 



Mp) 



^Z.Z,Vip,k)^l^^Aik) (4) 



(27r)3 



2Ek 



which is obtained as the quasiparticle limit of the gap 
equation with dressed propagators (|^). The superfluid 
energy gap is 



where 



and 



A{p) = ZpA{p) 

I _ dj:{p,ujp) 

Ep = ^ {ujp - Ai)2 + A{py 



(5) 



(6) 



(7) 



The usual BCS equation is obtained by putting Zp = 1 
in the above equations. In that case, the only influence 
of the medium comes through the modification of the 
dispersion relation ujp. 

Dispersive self-energy corrections to the gap equation 
are twofold. First, the superfluid energy gap is multiplied 
by the quasiparticle strength Zp^ < 1 (||), and second, 
the interaction strength in the gap equation is reduced 
by a factor ZpZ^ (^). In the case of symmetric nuclear 
matter at normal density the reduction of the superfluid 
gap due to self-energy effects is of about one order of 
magnitude J 2l[| . 

we present the results of an analogous cal- 
the pairing in neutron matter for two 
interaction potentials |2^. It is known that for the 
^5*0 gap different interaction potentials lead to similar 
superfluid gaps Q| . We recall that the spectral functions 
and the single-particle energies used in the gap equation 
are obtained from a separate self-consistent calculation 
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FIG. 2: The ^So superfluid energy gap in neutron matter, 
for the CD- Bonn interaction in the BCS approximation (soHd 
Une), from the gap equation with dressed proptisators (^) 
(dashed line), and from the effective gap equation (dashed 
line with stars); and for the Argonne V18 interaction in the 
BCS approximation (dashed-dotted line) , from the gap equa- 
tion with dressed propagators (^) (dotted line), and from the 
effective gap equation (W) (dotted line with stars). 



in the normal phase. We find the Zp factor at the Fermi 
momentum ~ 0.9 for the range of densities studied. It 
is close to the one obtained from the BHF calculation 
without rearrangement terms pC| |. The rearrangement 
terms in the BHF approach reduce this value further 
(the T-matrix result includes the rearrangement contri- 
bution to the self-energy). We find a smaller Z factor 
than the authors of reference The obtained values 
of Zp translate directly into the resulting superfluid gap, 
which we find in between the results of Refs. and 
pO| . The average reduction factor is 0.7 in the region of 
the maximal gap. The effective gap equation (||) is an 
excellent approximation to the solution of the gap equa- 
tion with full spectral functions (H). It is due to the fact 



that in this region of densities, the gap equation kernel is 
dominated by contributions close to the Fermi momen- 
tum, where the quasiparticle limit applies. We could not 
perform the self-consistent T-matrix calculation at den- 
sities below O.OlGfm^'^; the expected disappearance of 
self-energy effects in the gap equation at low densities 
could not be explicitly demonstrated. This is in con- 
trast to the results for the energy per particle where the 
BHF approach and the self-consistent T-matrix calcula- 
tion merge at low density ||2^ (Fig ^ . 

In the gap closure region the ratio of the BCS gap to 
the gap obtained with dressed propagators is large. Such 
a strong reduction of the gap was found in symmetric 
nuclear matter at normal density pTf , which is close to 
gap closure for the '^S'l Di gap. Also the effective 
gap equation (^) and the gap equation with full spectral 
functions (^ give markedly different results. 

In summary, we study pure neutron matter using the 
in-medium T-matrix approximation. We present a com- 
parison of the energy per nucleon for neutron matter in 
the BHF and T-matrix approaches. We find similar bind- 
ing energies at low densities, but a stiffer equation of state 
for the T-matrix calculation at higher densities. This is 
an explicit confirmation of the expected similarities and 
discrepancies between the two approaches. 

We use self-consistent spectral functions in neutron 
matter to calculate the superfluid gap from the gap equa- 
tion with dressed propagators. This allows us to test the 
accuracy of the effective gap equation (^), which incor- 
porates to some extend self-energy effects. Using full 
spectral functions, we find, in the maximal gap region, a 
reduction of the gap of about 30% and a very good agree- 
ment with the effective gap equation. In the gap closure 
region, we find a stronger reduction. In this region the 
effective gap equation (^) is a less reliable estimate of self- 
energy effects. It should be remembered that also other 
in medium effects modify the superfiuid gap in neutron 
26l EtI Esl. 
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